Robust creation of arbitrary-sized Dicke states of trapped ions by global addressing 
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We propose a novel technique for the creation of maximally entangled symmetric Dicke states 
in an ion trap using adiabatic passage, which requires only a pair of chirped pulses from a single 
laser and is applicable to any number of ions and excitations. By utilising a particular factorisation 
of the Hilbert space for multi-level ladders we show that the problem can be reduced to 'bow-tie' 
configuration energy-level crossings. This technique is naturally robust against fluctuations in the 
laser intensity and the chirp rate. Even when realistic heating rates are considered, we estimate 
that the overall fidelity should remain high (e.g. 98% for a 10- ion Dicke state), which represents a 
significant improvement over traditional approaches. 
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Introduction. Genuine multipartite entangled states 
involving many excitations are centrally important for 
quantum information science. They are essential for 
quantum communication across more than one chan- 
nel, provide a universal resource for quantum compu- 
tation and can be used to improve the sensitivity of 
atomic clocks. Furthermore, developing a clear under- 
standing of entanglement is fundamental to our descrip- 
tion of microscopic systems. However, both the produc- 
tion and characterisation of large entangled states remain 
formidable challenges. A particularly important class of 
multi-particle entangled states are the Dicke-symmetric 
states of N particles and m excitations: 
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where \Pk] denotes the set of all distinct combinations 
of ions [l|] and = Nl/ [m\(N — m)!]. These states are 
robust against decohcrcncc, particle loss and measure- 
ment. Also, the entanglement that they contain cannot 
be destroyed by local operations Q . 

Currently, the maturest experimental approach to the 
creation of large entangled states is to couple trapped 
ions by allowing them to interact with a common vibra- 
tional mode. This can be achieved either with resonant 
laser pulses of precise area 0, H| or off-resonant geometric 
gates 0, @], and Dicke states containing a single quan- 
tum of energy (the so-called VK-states) have been fabri- 
cated in chains of up to 8 ions using the first technique 
[J. However, addressing the ions one or two at a time 
requires a costly increase in the number of preparation 
steps as the complexity of the final state grows, and for 
this reason, Dicke-symmetric states with more than one 
quantum of energy represent a significant theoretical and 
experimental challenge. The creation of Dicke-states has 
been proposed using simultaneous addressing for a cou- 
pled spin-chain, by using many precisely-timed pulses to 
navigate through a network of energy-level crossings Q. 
In this Letter, we propose an extremely simple method 



for producing Wj^ ) states robustly, for arbitrary m and 
N using a two-step adiabatic-passage technique. 

Level-scheme. We consider N ions cooled and 
trapped in a linear array. Each ion has two relevant in- 
ternal states, |0) and |1), and interacts with a laser pulse 
that is tuned close to (but not exactly at) the frequency 
of the first red sideband of the lowest vibrational mode. 
We assume that the phonon spectrum can be resolved 
sufficiently well that only the levels of the centre of mass 
mode are excited by this interaction and that other vibra- 
tional modes can safely be neglected. In the Lamb-Dickc 
limit, the interaction Hamiltonian is (h = 1) 
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where fi^ is the Rabi-frequency for the j ion, r\ is 
the single-ion Lamb-Dicke parameter, Aj is the j th ion- 
laser detuning, and N the number of ions. Equa- 
tion ([2|) holds in the limit that the trap frequency, 
i/ 2 > {2m ] {t)ri) 2 /N @, %. Since the total number of 

excitations, N = afa + Ylj=i |l)j(l|j> i s conserved, the 
dynamics within the subspace corresponding to a given 
number of quanta is closed. It is sufficient to consider 
the case of m-excitations shared collectively between the 
TV ions and the bus mode 0]. This subspace splits nicely 
into manifolds corresponding to/z = 0,l,...,m phonons 
in the bus-mode, each of which is (C ? J[_ At )-fold degener- 
ate. The relevant level scheme is sketched in figure [U 



Morris-Shore factorisation. Even though the total 
number of energy quanta is preserved by the Hamilto- 
nian, ([2]) , the relevant Hilbert space contains J2™=o ^ 
states, which grows very rapidly as either the number 
of ions or the number of excitations is increased. This 
means that an exact analytic treatment of this sys- 
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FIG. 1: Level-scheme for N ions coupled to their common 
centre-of-mass mode and sharing m quanta. (The detuning is 
not shown in this figure.) The number of quanta is conserved 
by ([2|) and so dynamics within the subspace containing m ex- 
citations is closed. This subspace consists of (m + 1) separate 
manifolds, each coupled only to its neighbours. 
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FIG. 2: Illustration of the Morris-Shore factorisation for 
N — 4 and m = 2. (a) Level-scheme for 4 ions and 2 exci- 
tations. Different colour lines represent the coupling between 
the bus mode and the four different ions, (b) The same sys- 
tem in the Morris-Shore basis consists of a single (m + l)-level 
subsystem, some m-level subsystems, some (m — l)-level sub- 
systems etc, all the way down to a collection of dark states. By 
ensuring that the required initial and final states are confined 
to the (m + l)-level subspace, the dynamics can be simplified 
considerably, since the coupling between different subsystems 



tcm rapidly becomes unfeasible using conventional ap- 
proaches. From a practical perspective, the number of 
control steps required to navigate between different states 
of this large Hilbert space soon becomes prohibitively 
high, since the overall fidelity is degraded with each con- 
trol step. Fortunately, a very convenient transformation 
exists which greatly simplifies the problem [H, [nj. By 
making a suitable change of basis, it is possible to de- 
couple the relevant state space into a series of smaller sub- 
spaces. These subspaces are not connected by the Hamil- 
tonian and so the dynamics within them remains closed. 
This procedure is known as the Morris-Shore factorisa- 
tion, and is sketched schematically in figure [2] Where it 
exists, this factorisation can provide a remarkable simpli- 
fication of the dynamics of many-level systems. However, 
in the general case, its existence is not guaranteed, since 



the couplings between the neighbouring manifolds arc re- 
quired to satisfy a specific commutation relation. 

Below, we shall show that for the present system the 
Morris-Shore transformation does exist, and furthermore 
it is possible to confine the dynamics to a single ladder 
involving (m + l)-levcls. By inducing a suitable level 
crossing within this subspace it is possible to transfer 
population adiabatically between the lowest and highest 
states of this (to + l)-level chain. In this way, Dickc- 
symmctric states may be robustly created in a single step 
for arbitrary N and to. The laser pulses we shall choose 
involve the same Rabi frequency Sl'(t) and detuning A(t) 
for all ions, and for ease of notation, we define the effec- 
tive Rabi-frequency f2(£) = r]£lj(t)/y/N. Experimentally 
this is an appealing choice, since this interaction Hamil- 
tonian may be realised by shining a single laser pulse onto 
all of the N ions. We proceed by changing the interaction 
picture to one in which the Hamiltonian can be written 
in block-matrix form: 
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The detuning-blocks A M (i) are proportional to the 
identity operator of dimension C^_ M (which we shall la- 
bel with equal elements (2fi — to) A(t)/2. The block 
represents the coupling between the manifolds 
with [i and (//— 1) phonons. The multi-level Morris-Shore 
transformation exists if and only if the commutation re- 

lation [v^-xWl^Wrf^J^V^it)] = holds 
for all the intermediate manifolds, (i.e. \i = 1, . . . , to — 1) 
i. 

Now, since Hj(t) is symmetric with respect to in- 
terchange of any two ions, the size of the non-zero 
elements in any one coupling block are the same. 
The magnitude of the non-zero elements in V^ iAI _i(t) 
is ft(t)y/JI/2. Therefore, we define scaled coupling 
blocks, V^-i = 2V M , M _i(£)/r2(i) v //I, which are all time- 
independent. The elements of these scaled coupling 
blocks can only be or 1. The above commutation rela- 
tion can now be re- written as follows: 
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Properties of the coupling blocks. There are four rel- 
evant groups of coupling, sketched in figure Q] (a) Go- 
ing from the (/! — l)-phonon to the /i-phonon manifold 
involves converting an ionic excitation into a phonon. 
Since each state in the (/i — l)-phonon manifold contains 
(to — /i+ 1) excited ions, each upper state is connected to 
this many lower states, (b) Coupling from the /z-phonon 
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to the (fx — l)-phonon manifold involves converting a 
phonon into an ionic excitation. Since each state in the 
jU-phonon manifold contains (N — to + /it) unexcited ions, 
each lower state is connected to (N — m + fi) upper states. 
Similarly, there are (N — m + /j> + 1) couplings in group 
(c) and (to — /i) in group (d). The products V^-i „_! 
can now be computed. The diagonal elements are 
(N — to + fx) , while the values of the off-diagonal ele- 
ments of Vp^-iVp „_ x arc simply the number of common 
states in the (/i — l)-phonon manifold coupled to both of 
the corresponding states in the /x-phonon manifold. (This 
is represented schematically by the red/yellow dashed 
lines in figure [TJ) These values can only ever be one (and 
there are (to — n)(N — m + fi) of these per row) or zero. 
For example, the states | . . . 10 . . and | . . . 01 . . .)\fj) 
are both coupled to the state | . . . 11 . . — 1), but there 
are no other states to which they can both be coupled 
by Vfi^-i. A similar analysis can be used to calculate 
the matrix Vj +1 „V^ + i jM . Once again, the values of these 
off diagonal elements arc one if these two states differ by 
the position of a single ionic excitation (blue/pink dashed 
lines in figure [T]) and zero otherwise. These off-diagonal 
elements are the same as in V^—iVl „_ 1; while the di- 
agonal elements are (to — /i) , giving: 



V^+i >M + (iV-2m + 2/i)V (4) 



Clearly, the matrices V^-iVl „_! an d 
commute, and so the generalised Morris-Shore transfor- 
mation always exists for to excitations shared collectively 
between N ions and the phonon mode. 

Morris-Shore states. Since there is only one TO-phonon 
state, there is only one (m + l)-level ladder in the Morris- 
Shore factorisation, and as we are interested in adia- 
batically coupling the lowest and highest states of this 
chain, only these two states are actually relevant. How- 
ever, in order to be sure which of the Morris-Shore eigen- 
states in the zero-phonon manifold belongs at the top of 
the ladder, it is necessary to calculate all (m + 1) lev- 
els. These are uniquely defined by the condition that for 
each pL = 1, . . . ,to, the eigenstates of V^.^-iV^ fl _ 1 and 

u-i^u.ju-i (in neighbouring manifolds) have the same 
eigenvalue 0] . After noting that the sum of the elements 
in any row of cither of these two matrices is equal to 
(to — fx + 1)(N — to + ju) , we conclude that the completely 
symmetric states, Wj^_„) \(J>) (with fj, = 0, . . . , to) are 
the eigenstates of the (to + l)-level chain in the Morris- 
Shore basis. This subspace is de-coupled from all other 
states during the dynamics and, conveniently, the state 
that we would like to create, i.e. jW^), is at the top of 
this ladder. 

We note that the only states involved in the dynamics 
are eigenstates of the combined ionic pseudospin, J z = 
\ (|l)i(l|j — 1 0) j (0 1 j ) - This is a necessary feature, 

which is made clear by writing the Hamiltonian in the 



form Hi(t) oc ( J+ae^ + h.c.V with J+ = J^f = i |l)i(0|j 
Below, we show how it is possible to use adiabatic 
techniques to steer the system into this target state. 

Adiabatic creation of Dicke- symmetric states. 
From now on, we shall only be concerned with the 
(m + l)-level ladder in the Morris-Shore factorisation. 
The Rabi-frequency fi(i) and detuning A(t) are se- 
lected in order to satisfy the objective of a smooth 
transfer of population between the two ends of this 
chain. Since the detunings, A M (i), are all proportional 
to A(t), chirping the laser across resonance will lead 
to a 'bow-tic' model, where all of the energies cross at 
the same instant of time. The relevant level scheme 
is shown in figure [3) Analytic bow-tie solutions have 
been found for three [ll[ and to states [l2| . The cou- 
plings between each pair of successive Morris-Shore states 
in the Dicke ladder with \i and fj, — 1 phonons are 
given by X^-i{t) = ^f-^fi(m - (J, + 1)(N - m + n). 
We choose to use a 'complex-sech' pulse, i.e. 
Q.(t) = fl a scch(t/T) and A(t) = A tanh(t/T). Com- 
pared to other types of level-crossing (e.g. Landau- 
Zener), these couplings result in a relatively fast adia- 
batic passage, with a very appropriate excitation profile; 
the excitation probability very close to unity in a certain 
range (controlled by Ao) around resonance, and vanish- 
ing exponentially rapidly outside this range [l3| . This is 
particularly relevant here, since it is important that our 
laser pulse does not excite unwanted vibrational modes. 
A transfer efficiency of 1 — e between the zero-phonon 
and TO-phonon states (with e <^ 1) requires the coupling 
strength and detuning to satisfy 
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FIG. 3: The m-phonon state and the Wj^} state can be 
smoothly connected adiabatically using a level-crossing, pro- 
viding that the couplings are chosen to satisfy the adiabaticity 
condition, ([5]). 

The Dicke-symmetric state, |Wj^), may be created as 
follows: (i) The system is prepared in its ground state 
and to of the ions (arbitrarily chosen) are promoted from 
|0) to This initialisation step may be performed ro- 
bustly, using adiabatic techniques [14|. (ii) A chirped 
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laser pulse is focused onto the m excited ions, creating an 
effective m-ion system containing m quanta. The other 
(N — m) ions remain in their ground states and de-couple 
from the dynamics. Providing the adiabaticity condition, 
([5]), is satisfied, the system will be transferred smoothly 
into the state |0 . . . 0) \m) . Thus our scheme also provides 
a method for adiabatically creating number states of the 
phonon field, which will be described in detail elsewhere 
1 1 51 ] . (Alternatively, this preparation may be achieved by 
performing stages (i) and (ii) on a single ion, m times.) 
(iii) Finally, the laser beam is broadened in space to cou- 
ple equally to all N ions. By applying a second chirped 
pulse, the system can be transferred smoothly into the 
target state, |W^}. 
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FIG. 4: The first three frames show the evolution of the popu- 
lations of all 22 states for the creation of a \ W%) state. In the 
frames labelled 1-phonon manifold and 0-phonon manifold, 
several indistinguishable curves are overlaid, since certain ions 
are equivalent. At time t ~ 8T, the state |Wi) is partially 
populated in the 1-phonon manifold. The bottom-right frame, 
labelled 'fidelity', shows the time-dependent overlap with the 
required final state (final fidelity = 99.996%). Parameters are 
QoT = 10; A T = 6. 

Figure [4] shows the result of a numerical simula- 
tion for the case of N = 6, m = 2. During the pe- 
riod — 10T < t < 0, population is transferred from state 
|1, 1,0,0,0,0)|0) to state |0,0,0,0,0,0)|2) via the 1- 
phonon manifold by addressing just 2 ions. After t = 0, 
the laser beam is spatially-broadened to couple equally 
to all 6 ions and the system is transferred back down 
into the 0-phonon manifold and into the symmetric state 
| )• Since our approach uses adiabatic passage to nav- 
igate between states, it is robust against fluctuations in 
the intensity and the timing of the laser pulses. Fur- 
thermore, even when intensity of the pulse was allowed 
to fluctuate by 10% across the chain, the overall fidelity 
remained above 99.3% for the parameters in figured] 

Heating effects. fl(t) is bounded from above by ap- 
proximately i/ tr /10, in order to prevent extraneous ex- 
citations in vibrational modes other than the centre of 
mass mode @. Together with the adiabaticity require- 



ment, ([5]), this confines the timescale, T for the chirped 
pulses to be at least 25/xs for a trap-frequency of 4MHz 
(as in [l6]). A conservative estimate for the time required 
to implement the entire technique is therefore 400/is (i.e . 
16T). Taking a phonon heating-rate of 5s _1 per ion [16[ 
we therefore estimate that for a 10-ion Dicke-state, the 
total amount of heating during the whole procedure could 
be on the level of 0.02 phonons; the infidelity introduced 
by bus- mode heating is thus of order 2%, which is sig- 
nificantly lower than could be expected from traditional 
approaches involving many sequential operations. More- 
over, the effects of heating are largely insensitive to the 
complexity the target state. 

Conclusions. We have shown that maximally- 
entangled Dicke states, |W)^), of arbitrary size may be 
generated in an ion-trap using a pair of extremely sim- 
ple laser pulses. Starting with m ions in state |1) and 
all other ions in state |0), the first chirped pulse simulta- 
neously addresses these m ions to create adiabatically a 
phonon number state with zero ionic excitations and m 
phonons. A second chirped pulse addresses all N ions to 
create the state |W^). This extremely simple scenario 
derives from three theoretical findings: (i) the multi-state 
Morris-Shore factorisation de-couples the vast Hilbcrt 
space of the N ionic qubits into independent smaller 
chainwisc linkages; (ii) the longest of these, an (m + 1)- 
dimcnsional chain, happens to contain all the Dicke states 
from \Wi) to IW 7 ^), along with the phonon-number 
state, \Wq ); (iii) using a frequency-chirped pulse cre- 
ates a bow-tie linkage pattern, wherein the state W(f ) 
connects adiabatically to the \W^) state. The proposed 
technique is adiabatic in nature and hence it is robust 
against intensity and frequency imperfections; moreover, 
because only two interaction steps are required, heating 
effects are reduced. We thank Peter Ivanov and Peter 
Blythc for helpful discussions. This work has been sup- 
ported by the EU ToK project CAMEL and the EU RTN 
project EMALI. 
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